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Abstract 

Synapsin 1 is phosphoprotein consist onto the surface of synaptic vesicles. It 
is involved in development of neurons. SYN1 play an important role in 
axonogenesis and synaptogenesis neurotransmitter release and connection of 
vesicles to the cytoskeleton. Mutations in SYN1 can cause schizophrenia. The 
alterations in SYNI lead to a _ neurodegenerative disorder known 
schizophrenia affects 1% of the population and is classified by lack of social 
behavior and various behavioral abnormalities. In silico approaches were 
employed to predict 3D structures and protein-protein docking of SYNI1. 
Comparative modeling and threading computational approaches were applied 
to predict 3D structure of SYN1. STRING database calculated 1AUV as the 
interacting partner of SYN1. SYN1 protein-protein molecular docking study 
was done against 1AUV to identify the interacting residues. The observed 
interacting residues may have effective response against schizophrenia. The 
generated findings of this effort help to design the effective therapeutic 
targets against SYN1 by targeting schizophrenia. 
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Introduction 


Synapsin 1 (SYN1) is a member of a phosphoprotein 
and also has other members Synapsin 2 and Synapsin 
3 [1, 2]. The synapsins (SYNs) exist in the central 
and peripheral nervous system that consist onto the 
top of synaptic vesicles [2, 3]. The family members 
have been identified by regular domain of protein, 
interlaced in neuronal expansion, involved in 
synaptogenesis and the transition of neurotransmitter 
break [4]. SYNs manage synaptic vesicles and 
elaborate to the management of synaptic vesicles [5, 
6]. 

The phosphoprotein SYN1 involves in axonogenesis 
and synaptogenesis regulation as well as in synaptic 
transmission and plasticity of mature neurons [5-7]. 
In vertebrate species, SYN1 plays an important role 
in the breakage of neural transmitter, axonogenesis 
and synaptogenesis demonstrated by mouse 
introducing with an epileptic constitution [5-8]. In 
nerve ending, SYNI participate in connecting the 
cavity to the cytoskeleton and to one another [9, 10]. 

SYNI involved in the development of neurons and 
also in the formation of synaptic junction between 
the neurons [11, 12]. SYN1 has two isoforms as 
isoform IA and IB with different lengths. The length 
of isoform JA is 705 and IB is 669 residues. During 
brain development, the expression of SYNI 
moderately increased and reaches at the highest level 
in adult neural tissues [13]. SYN1 is phosphor- 
protein plays multiple roles in the synaptic 
transmission [14-16]. The SYN1 is phylogenetically 
conserved and has been mapped on chromosome X 
in human, mouse and other mammals [17, 18]. 

Non sense and missense mutations in human SYN1 
cause epilepsy and autism. Autism disorder occurred 
in the presence and absence of epilepsy [19, 20]. The 
modifications in SYNI are due to mutations, that 
causing imperfection in synaptic vesicles and nerve 
terminal functions [21]. According to its natural 
purpose, SYNI is present in brain and neuron 
represented by the supporter domain of SYNI [22]. 
The SYNI act as a substance for various protein 
phosphokinases and phosphorylation and also 
perform key role in the management of SYNI in 
nerve ending [23-26]. 

Garcia et al., explained that only variations causing a 
imperfection in SV trade and neural transmitter 
breakage in humans [27], announced a W356X 
modification in SYNI in a family separating a 
jumble of neural expansion defects along with 
epilepsy. SYNs are a family member of neuron 
specific SV phosphoproteins interlaced in synaptic 
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transmission and softness. SYNs manage SV trade in 
the middle of reserve pool (RP) and readily 
releasable pool (RRP). SYNs also involved in 
neuronal evolution, synaptogenesis and preservation 
of adult synapses [14, 28]. It has been recognized 
that variations (one non-sense and three mis-sense) 
within SYN1 connected with Epilepsy and Autism 
Spectrum Disorders. Three out of four variations 
were gathered in SYN1I D-domain that were related 
to effective flaws in nerve ending function [29]. 
Schizophrenia is a severe debilitating psychiatric 
mental disorder with heritability about 80% 
unorganized speeches, abnormal behave and illusion 
of imagination. Schizophrenia is characterized by 
abnormal behavior including obsession, dampening 
of emotions and self-deceptions. Approximately 1% 
of the population effected from schizophrenia in late 
teen to early twenty period of lifetime. Many 
analyses provided helping information to the plan 
that malfunction of SYN are associated to the 
schizophrenia onset. Over years some data analyses 
assembled characterized expression level of different 
SYNs on post-mortem brain sample of human. 
Display general decreasing of the size of SYNs 
protein and various expressions of SYNs 
modifications in brains obtain from schizophrenia. 
Bioinformatics is an interdisciplinary field of science 
helps to solve numerous biological problems by 
utilizing molecular biology, computational science 
and statistics. Computer aided drug design has 
progressive success in computational research 
methodologies to solve biological problems [30-32] 
and computationally designed molecules against 
cancer [33, 34] and neurological disorders [35, 36]. 
The present study has structural modeling and 
identification of drug target sites of SYNI against 
schizophrenia. The NMR and X-ray crystallographic 
structure of SYN1 has not been reported yet. 


Material and Methods 


Structural bioinformatics involved in the prediction 
of protein structure. In case of neurological diseases, 
it has constructive methods to sketch functional 
composite. SYNs involved in development of 
neurons. The expression of SYNs_ in brain 
development cause nervous disorders. 

SYNI protein sequence was retrieved from Uniprot 
KB having the accession number P17600. In this 
work, structure prediction and protein-protein 
docking analyses were executed. 

The amino acid sequence of SYN1 was subjected to 
BLASTp for the selection of suitable templates 
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against Protein Data Bank (PDB) [37]. MODELLER 
9.20 was utilized to predict the 3D structure of 
SYN1. 3D structures of SYN1 were predicted from 
MODELLER 9.20 [38]. For structure projection, 
SWISS MODEL [39], I-TASSER, MOD-WEB, 
Phyre2, HHpred, PSIPRED, CPHModel, M4t [40], 
Robetta [41], SPARKS-X, Raptor-X were utilized. 
UCSF Chimera 1.13 was used for the visualization of 
SYN1 3D structures. Anolea [42], ERRAT [43], 
Rampage and Verify 3D [44] were used structure 
evaluation. The graph and excel sheet described 
overall assessment of the protein structures. For the 
analyses of overall quality factor, Errat software was 
employed. Anolea evaluation tool was used for the 
calculation of Z-score of all predicted structures. 
Rampage tool was used to calculate the favored 
region, allowed region and outliers. The energy 
minimization was done by UCSF Chimera 1.13.1 
[45]. 

The geometry expansion and energy minimization of 
the molecules were executed by Chem3D Ultra. 
PatchDock, FireDock, and PyRx were utilized for 
comparative molecular docking. Ligplot and UCSF 
Chimera 1.13.1 were used for visualization. 
STITCH4 [46] (Search Tool for InTeracting 
CHemicals) and STRING [47] (Search Tool for the 
Retrieval of INteracting Genes/Proteins) were used 
to examine the interacting partner of SYNI. 
PatchDock was used for verifying as well as 


Table 1: The selected templates through BLASTp. 
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acceptance of the generated protein-protein results 
[48, 49]. 


Result and Discussion 


Extensive literature review and biological databases 
showed interesting information about the sequence of 
SYNI. Ala, Pro and Ser were observed polar 
residues perform significant functions including 
stability, folding and function of SYN1 (Fig. 1). 


Structure Prediction 


Comparative modeling and threading were used to 
predict 3D structures of SYN1. The sequences of 
SYNI were subjected to BLASTp against PDB 
database used to search suitable templates (Table 1). 
The five suitable templates with properties of query 
coverage, maximum identity, E values and total score 
were selected for homology modeling. The five 
selected templates were used to generate 50 
structures in which each template generated 10 
structures. The templates were used to predict 3D 
structures of SYN1. The highest query coverage of 
suitable templates was 59% and leftover templates 
have query coverage less than 42%. All the generated 
models were evaluated for favored region, allowed 
region, Outliers, overall quality factor and Z-Score. 
By comparison of generated models, a comparative 
graph was generated (Fig. 2). 


Description Max Score Total Score Query Cover E-Value Per. Ident Accession 
Chain A, Rat synapsin I 864 864 59% 0.0 97.14% 1PK8 
Chain A, Synapsin Ia 670 670 44% 0.0 98.39% 1AUX 
Chain A, Synapsin Ia 637 637 44% 0.0 94.86% 1AUV 
Chain A, crystal Structure Analysis of the 543 543 43% 0.0 77.92% LI7L-A 
complex of the C domain of Synapsin | 

Chain A, SYnapsin-3 498 498 42% le-172 70.72% 2POA 
ERRAT was utilized to evaluate quality factor of Protein-protein dockin g studies 

SYNI. The energy minimization on model forecast 

por OE SINI WaS APP lied to refine the The protein-protein docking analyses were 
stereochemistry. The suitable structure was 


examined. The chosen structures after assessment 
restriction were subjected to UCSF Chimera 1.10 for 
minimization at 1000 steepest and conjugate slops 
runs (Fig. 3). String was applied to predict 
interacting part of SYN1. The structure has highest 
identity subjected with UCSF Chimera 1.10 and 
removed extra residues. 


performed through PatchDock. The receptor ligand 
(SYN1-1AUV) complex was examined onto the base 
of atomic contact energy and top ranked 10 
complexes. Additionally, analyses were evaluated 
through FireDock and the complexes were examined 
on to the base of minimum binding global energy 
(Table 2) (Fig. 4). 
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Fig. 1: Pie chart representation of composition amino acids and calculated percentage values of SYN1. 
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Fig. 2: Graph shows the quality factor (orange brown), allowed region (yellow), favored region (blue) and outlier 
region (red) “From all structures, most reliable structure was selected from graph. 
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Fig. 3: The protein 3D structure of SYN1. 
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Fig. 4: The protein-protein docking structure analyses. 
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Table 2: The interacting residues of target protein and interacting protein 
Target protein Target Protein Residues Interacting Protein Interacting Protein Residues 


SYN I GLU 289, GLY 287, PHE 492, PRO 262, 1AUV ARG 381, ASP 379, VAL 331, MET 343, ALA 342, 
ALA 491, THR 579, ALA 580, SER 577 VAL 280, MET 270ILE 291, VAL 294 
P [8] Humeau Y, Candiani S, Ghirardi M, Poulain B, 
Conclusion Montarolo P. Functional roles of synapsin: lessons 


Schizophrenia is considered as a complex disorder 
affects almost 1% of the total population. 
Computational techniques were utilized based on 
structural modeling docking analyses. The 3D 
structure prediction proposed the reliability of the 
structure. The docking analyses satisfied the least 
binding global energy and elucidating the target 
binding sites of SYN1 against schizophrenia. 


Acknowledgments 


The authors are thankful to Mr. Jonathan Javid, for 
their contribution and guidance for conducting this 
research. 


Conflict of interest 


The authors declare no conflicts of interest. 


References 


[1] Patzke C, Brockmann MM, Dai J, Gan KJ, Grauel MK, 
Fenske P, et al. Neuromodulator signaling 
bidirectionally controls vesicle numbers in human 
synapses. Cell. 2019;179:498-513. e22. 

[2] Hedegaard C, Kjaer-Sorensen K, Madsen LB, 
Henriksen C, Momeni J, Bendixen C, et al. Porcine 
synapsin 1: SYNI gene analysis and functional 
characterization of the promoter. FEBS open bio. 
2013;3:411-20. 

[3] Fornasiero EF, Bonanomi D, Benfenati F, Valtorta F. 
The role of synapsins in neuronal development. 
Cellular and molecular life sciences. 2010;67:1383-96. 

[4] Farisello P, Boido D, Nieus T, Medrihan L, Cesca F, 
Valtorta F, et al. Synaptic and extrasynaptic origin of 
the excītation/inhibition imbalance in the 
hippocampus of synapsin I/II/HI knockout mice. 
Cerebral cortex. 2013;23:581-93. 

[5] Shupliakov O, Haucke V, Pechstein A. How synapsin I 
may cluster synaptic vesicles. Seminars in cell & 
developmental biology: Elsevier; 2011. p. 393-9. 

[6] Nemani VM, Lu W, Berge V, Nakamura K, Onoa B, 
Lee MK, et al. Increased expression of a-synuclein 
reduces neurotransmitter release by inhibiting synaptic 
vesicle reclustering after endocytosis. Neuron. 
2010;65:66-79. 

[7] Hart A. Role of synapsin in long-term synaptic 
facilitation in Aplysia. 2011. 


from invertebrates. Seminars in cell & developmental 
biology: Elsevier; 2011. p. 425-33. 

[9] Tyebji S, Hannan AJ. Synaptopathic mechanisms of 
neurodegeneration and dementia: Insights from 
Huntington’s disease. Progress in neurobiology. 
2017;153:18-45. 

[10] Michael R. Analysis of Synapsin I and II mRNAs and 
Proteins During Murine Development. Mol Biol. 
2017;6:2. 

[11] Li L, Xiong W-C, Mei L. Neuromuscular junction 
formation, aging, and disorders. Annual review of 
physiology. 2018;80:159-88. 

[12] Martinez G, Khatiwada S, Costa-Mattiol M, Hetz C. 
ER proteostasis control of neuronal physiology and 
synaptic function. Trends in neurosciences. 
2018;41:610-24. 

[13] Bogen I, Jensen V, Hvalby O, Walaas S. Synapsin- 
dependent development of glutamatergic synaptic 
vesicles and presynaptic plasticity in postnatal mouse 
brain. Neuroscience. 2009;158:231-41. 

[14] Cesca F, Baldelli P, Valtorta F, Benfenati F. The 
synapsins: key actors of synapse function and 
plasticity. Progress in neurobiology. 2010;91:313-48. 

[15] A. Sehgal SAFK, Seemab. In Silico Analyses, 
Bioequivalence and Disposition Kinetics of 
Allopurinol in Healthy Male Subjects. Drug Delivery 
Letters. 2016;6:113-21. 

[16] By Sheikh Arslan Sehgal AHM, Rana Adnan Tahir, 
Asif Mir. Quick Guideline for Computational Drug 
Design. Bentham Science Publisher. 2018. 

[17] Kügler S, Kilic E, Bähr M. Human synapsin 1 gene 
promoter confers highly neuron-specific long-term 
transgene expression from an adenoviral vector in the 
adult rat brain depending on the transduced area. Gene 
therapy. 2003;10:337-47. 

[18] Tahir RA, Sehgal SA. Pharmacoinformatics and 
Molecular Docking Studies Reveal Potential Novel 
Compounds Against Schizophrenia by Target SYN II. 
Comb Chem High Throughput Screen. 2018;21:175- 
81. 

[19] Garcia C, Blair H, Seager M, Coulthard A, Tennant S, 
Buddles M, et al. Identification of a mutation in 
synapsin I, a synaptic vesicle protein, in a family with 
epilepsy. Journal of medical genetics. 2004;41:183-6. 

[20] Fassio A, Patry L, Congia S, Onofri F, Piton A, 
Gauthier J, et al. SYN1 loss-of-function mutations in 
autism and partial epilepsy cause impaired synaptic 
function. Human molecular genetics. 2011;20:2297- 
307. 

[21] Thiel G, Ekici M, Rössler OG. RE-1 silencing 
transcription factor (REST): a regulator of neuronal 
development and neuronal/endocrine function. Cell 
and tissue research. 2015;359:99-109. 


202 


[22] Schoch S, Cibelli G, Thiel G. Neuron-specific gene 
expression of synapsin I major role of a negative 
regulatory mechanism. Journal of Biological 
Chemistry. 1996;271:3317-23. 

[23] Kugler S, Meyn L, Holzmiiller H, Gerhardt E, 
Isenmann S, Schulz J, et al. Neuron-specific 
expression of therapeutic proteins: evaluation of 
different cellular promoters in recombinant adenoviral 
vectors. Molecular and Cellular Neuroscience. 
2001;17:78-96. 

[24] Scott BB, Lois C. Generation of tissue-specific 
transgenic birds with lentiviral vectors. Proceedings of 
the National Academy of Sciences. 2005;102:16443-7. 

[25] Hioki H, Kameda H, Nakamura H, Okunomiya T, 
Ohira K, Nakamura K, et al. Efficient gene 
transduction of neurons by lentivirus with enhanced 
neuron-specific promoters. Gene therapy. 
2007;14:872-82. 

[26] Sims K, Ahmed Z, Gonzalez AM, Read ML, Cooper- 
Charles L, Berry M, et al. Targeting adenoviral 
transgene expression to neurons. Molecular and 
Cellular Neuroscience. 2008;39:411-7. 

[27] Goodship J, Garcia C, Blair H, Seager M, Coulthard 
A, Tennent S, et al. Deficiency of synapsin I, a 
Synaptic vesicle protein, causes epilepsy. Journal of 
Medical Genetics: British Human Genetics 
Conference: Newcastle University; 2003. 

[28] Greco B, Manago F, Tucci V, Kao H-T, Valtorta F, 
Benfenati F. Autism-related behavioral abnormalities 
in synapsin knockout mice. Behavioural brain 
research. 2013;251:65-74. 

[29] Lo-Castro A, Curatolo P. Epilepsy associated with 
autism and attention deficit hyperactivity disorder: is 
there a genetic link? Brain and Development. 
2014;36:185-93. 

[30] Kanwal S, Jamil F, Ali A, Sehgal SA. Comparative 
Modeling, Molecular Docking, and Revealing of 
Potential Binding Pockets of RASSF2; a Candidate 
Cancer Gene. Interdiscip Sci. 2017;9:214-23. 

[31] Sehgal SA. Pharmacoinformatics, Adaptive Evolution, 
and Elucidation of Six Novel Compounds for 
Schizophrenia Treatment by Targeting DAOA (G72) 
Isoforms. Biomed Res Int. 2017;2017:5925714. 

[32] Sehgal SA. Pharmacoinformatics and molecular 


docking studies reveal potential novel Proline 
Dehydrogenase (PRODH) compounds for 
Schizophrenia inhibition. Medicinal Chemistry 


Research 2017;26:314-26. 

[33] Sehgal SA, Hammad MA, Tahir RA, Akram HN, 
Ahmad F. Current Therapeutic Molecules and Targets 
in Neurodegenerative Diseases Based on in silico 
Drug Design. Curr Neuropharmacol. 2018;16:649-63. 

[34] Sehgal SA, Khattak NA, Mir A. Structural, 
phylogenetic and docking studies of D-amino acid 
oxidase activator (DAOA), a candidate schizophrenia 
gene. Theor Biol Med Model. 2013;10:3. 


Biomedical Letters 2020; 6(2):197-203 


[35] Sehgal SA, Mannan S, Kanwal S, Naveed I, Mir A. 
Adaptive evolution and elucidating the potential 
inhibitor against schizophrenia to target DAOA (G72) 
isoforms. Drug Des Devel Ther. 2015;9:3471-80. 

[36] Sehgal SA, Mannan S, Ali S. Pharmacoinformatic and 
molecular docking studies reveal potential novel 
antidepressants against neurodegenerative disorders by 
targeting HSPB8. Drug Des Devel Ther. 
2016;10:1605-18. 

[37] Rose PW, Prlić A, Altunkaya A, Bi C, Bradley AR, 
Christie CH, et al. The RCSB protein data bank: 
integrative view of protein, gene and 3D structural 
information. Nucleic acids research. 2016:gkw1000. 

[38] Webb B, Sali A. Protein structure modeling with 
MODELLER. Protein structure prediction: Springer; 


2014. p. 1-15. 

[39] Arnold K, Bordoli L, Kopp J, Schwede T. The 
SWISS-MODEL workspace: a web-based 
environment for protein structure homology 


modelling. Bioinformatics. 2006;22:195-201. 

[40] Fernandez-Fuentes N, Madrid-Aliste CJ, Rai BK, 
Fajardo JE, Fiser A. M4T: a comparative protein 
structure modeling server. Nucleic acids research. 
2007;35:W363-W8. 

[41] Kim DE, Chivian D, Baker D. Protein structure 
prediction and analysis using the Robetta server. 
Nucleic acids research. 2004;32:W526-W31. 

[42] Melo F, Devos D, Depiereux E, Feytmans E. 
ANOLEA: a www server to assess protein structures. 
Ismb1997. p. 187-90. 

[43] Dym O, Eisenberg D, Yeates T. ERRAT. 2012. 

[44] Grotthuss Mv, Pas J, Wyrwicz L, Ginalski K, 
Rychlewski L. Application of 3D-Jury, GRDB, and 
Verify3D in fold recognition. Proteins: Structure, 
Function, and Bioinformatics. 2003;53:418-23. 

[45] Pettersen EF, Goddard TD, Huang CC, Couch GS, 
Greenblatt DM, Meng EC, et al. UCSF Chimera—a 
visualization system for exploratory research and 
analysis. Journal of computational chemistry. 
2004;25:1605-12. 

[46] Kuhn M, von Mering C, Campillos M, Jensen LJ, 
Bork P. STITCH: interaction networks of chemicals 
and proteins. Nucleic acids research. 2007;36:D684- 
D8. 

[47] Franceschini A, Szklarczyk D, Frankild S, Kuhn M, 
Simonovic M, Roth A, et al. STRING v9. 1: protein- 
protein interaction networks, with increased coverage 
and integration. Nucleic acids research. 
2012;41:D808-D15. 

[48] Tahir RA, Hassan F, Kareem A, Iftikhar U, Sehgal 
SA. Ligand-Based Pharmacophore Modeling and 
Virtual Screening to Discover Novel CYP1A1 
Inhibitors. Curr Top Med Chem. 2019;19:2782-94. 

[49] Tahir RA, Bashir A, Yousaf MN, Ahmed A, Dali Y, 
Khan S, et al. In Silico identification of angiotensin- 
converting enzyme inhibitory peptides from MRJP1. 
PLoS One. 2020;15:e0228265. 


203 


